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LONG-TERM  STABILITY  FOR  PARTICLE  ORBITS 

This  report  contains  the  lecture  delivered  by  Dr.  L.  J.  Laslett 
at  a  one-day  conference  on  The  Mathematical  and  Computational 
Aspects  of  Accelerator  Design,  and  a  discussion  by  Dr.  J.  Moser 
of  some  of  the  points  raised  by  Dr.  Laslett.   Two  appendices 
contain  some  further  comments  by  Dr.  Laslett  and  Dr.  E.  McMillan. 
In  addition  to  the  talks  reported  here,  the  Conference  also 
heard  a  report  by  Dr .  G.  Parzen  on  the  iterative  techniques 
used  to  calculate  the  eigenfrequencies  in  recent  cavities  in 
connection  with  linear  accelerator  design.   He  reported  on  the 
difficulties  caused  by  shapes  that  consist  of  two  nearly  dis- 
connected domains  of  unequal  size,  and  of  methods  for  over- 
coming these  difficulties.   Dr.  Symon  reported   on  non-linear 

resonance  phenomena  due  to  the  coupling  of  two  particles,  and 
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Dr.  E.  Courant  reported   on  the  effect  of  coupling  in  many 

particles.  Dr.  A.  Garren  described  some  sophisticated  features 
of  automation  in  codes  used  in  accelerator  design. 

We  thank  all  participants  of  the  Conference  for  their 
contribution,  and  hope  that  this  report  will  stimulate  mathemati- 
cians to  concern  themselves  with  aspects  of  long-range  stability. 

See  H.  Meier  and  K.  R.  Symon,  "Analytical  and  Computational 
Studies  on  the  Interaction  of  a  Sum  and  a  Difference  Resonance," 
Proc .  of  the  Intern.  Conf.  on  High-Energy  Accelerators  and 

Instrumentation  -  CERN  1959 
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See  E.  D.  Courant  and  A.  M.  Sessler,  "Transverse  Coherent 

Resistive  Instabilities  of  Azimuthally  Bunched  Beams  in 

Particle  Accelerators,"  The  Review  of  Scientific  Instruments, 

Vol.    37,    No.    11,    1579-1588,    Nov.,    1966. 
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LONG-TERM  STABILITY  FOR  PARTICLE  ORBITS 

■X- 

L.  Jackson  Laslett 


The  three  accelerator  papers  on  the  program  for  this 
afternoon  will  be  concerned  with  types  of  dynamical  problems 
that  arise  in  the  study  of  particle  accelerators,  although 
the  same  mathematical  problems  may  well  arise  also  in  other 
contexts.   The  accelerators  with  which  we  are  concerned 
presumably  will  be  of  the  "alternating-gradient"  type,  in 
which  --  for  very  good  reasons  --  the  functional  character 
of  the  focusing  force  experienced  by  an  individual  displaced 
particle  changes  periodically  as  the  particle  proceeds  on  its 
way.   Dr.  Symon  will  review  some  of  the  quite  striking  effects 
that  can  develop  from  " coupling  resonances,"  when  the  equations 
of  motion  are  non-linear,  and  Dr.  Courant  will  summarize  the 
complications  that  develop  with  beam  intensities  sufficiently 
great  that  inter-particle  forces  must  be  considered. 

My  own  paper  is  intended  to  report  results  on  some 

computer  experiments  that  I  had  hoped  might  cast  light  on 

the  question  of  long-term  stability  --  specifically,  in  my 

work,  for  motion  with  only  one  spatial  degree  of  freedom. 

For  most  accelerator  projects  this  question  may  not  be  one 

of  as  immediate  importance  as  the  questions  to  be  discussed 

by  my  colleagues  this  afternoon,  but  the  answers  could  be 

vital  for  the  proton  storage-ring  devices  now  being  considered 

*  Lawrence  Radiation  Laboratory,  University  of  California, 
Berkeley,  California 


or  under  construction.   The  general  problem  appears  to  be  a 
very  interesting  and  quite  difficult  one,  and  I  would  look 
forward  to  comments  from  our  mathematician  friends  concerning 
what  one  can  say  at  this  time  in  simple  terms  with  respect 
to  this  problem. 

Thus,  for  my  own  part,  I  would  like  this  afternoon  to 
call  to  your  attention  this  question  of  long-term  stability 
and  report  on   a  few  elementary  computational  experiments  that 
I  performed  last  summer  for  my  own  orientation  and  amusement. 
As  for  motivation,  the  first  Figure  indicates  some  of  the 
numbers  that  may  be  relevant  for  characterizing  modern 
accelerators  with  respect  to  the  interval  over  which  one  would 
wish  them  to  exhibit  stabij.ity  of  the  particle  motion.   We 

note  in  particular  that  a  representative  particle  may  be  called 
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upon  to  traverse  some  10   periods  of  the  alternating-gradient 

(A-G)  focusing  structure.   These  parameters  are  indeed 

essentially  the  same  as  the  number  of  oscillations  or  periods 

successfully  experienced  by  the  beams  of  existing  high-energy 

synchrotrons,  and  so  we  have  the  opportunity  of  concluding  from 

experience  that  treacherous  long-term  instabilities  need  not 

be  troublesome  in  a  conventional  A-G  design. 


PRESENT  AND  FUTURE  A-G  SYNCHROTRONS 

Oscillation  Wavelength: One  to  Several  Hundred  Meters 
Period  of  Structure:     10  to  100  m. 

Time  Interval:  ^   1  sec. 

No.  of  Periods  Traversed: 


(3  xl0°  m/sec) (1  sec)    , „7 
(10  to  100  m) 


Fig.  1 

It  would  not  "be  out  of  the  question,  moreover,  to  think  of 
computational  experiments  that  would  subject  particles  with 

a  limited  number  of  selected  initial  conditions  to  algebraic 
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transformations  that  could  simulate  passage  through  some  10 

periods  of  the  A-G  structure,  although  one  would  have  to  be 
attentive  to  the  possibility  of  obtaining  misleading  results 
as  a  consequence  of  round-off  or  truncation  errors  --  especially 
if  the  structure  of  the  problem  is  such  that  the  distinction 
between  the  initial  conditions  for  stable  vs.  unstable  motion 
is  very  fine  grained, 

A  storage  ring,  on  the  other  hand,  might  be  of  comparable 


dimensions  --  or  at  least  Involve  a  similar  number  of 
oscillations  per  second  and  number  of  periods  traversed  per 
second  --  but  would  be  intended  to  retain  particles  for  hours . 

Thus,  as  we  see  on  Fig.  2,  one's  interest  in  long-range  stability 

4 
becomes  extended  to  intervals  some  10   times  as  great  --  i.e., 

possibly  to  10   periods,  or  even  to  5  x  10   periods  of  the 

structure . 


PRESENT  AND  FUTURE  A-G  SYNCHROTRONS 

Oscillation  Wavelength: One  to  Several  Hundred  Meters 
Period  of  Structure:     10  to  100  m. 

Time  Interval:  ~  1  sec. 

No.  of  Periods  Traversed: 


Q 

(3  X  10  m/sec)  (1  sec)   ^  -^^7 
(10  to  100  m) 


STORAGE  RINGS 

k 

Time  Interval:  ~  3  hrs.  or  'v  10   sec 

No.  of  Periods  Traversed:  ^   10 


Fig.  2 


Here  past  experience  is  not  particularly  helpful,  since  the 
storage  rings  that  have  been  operated  to  date  have  been  electron 
rings  (or  electron-positron  rings)   in  which  radiation  damping 
can  helpfully  play  a  dominant  role  that  is  effectively  absent 
(by  a  relative  factor  of  ten  orders-of -magnitude  or  more)  for 
protons . 

Cosmologically,  the  universe  is  supposed  to  be  only 

Q 

some  5  x  10  years  old,  corresponding  to  no  more  than  that  of 
many  planetary  periods  at  the  present  rate,  so  astronomical 
evidence  may  provide  little  assurance  to  the  storage-ring 
designer  unless  we  are  prepared  to  be  relatively  restrained 
in  our  non-linearities,  and  unless  we  have  some  theoretical 
scaling  laws  that  then  could  be  applied  to  relate  the  long- 
term  stability  of  our  proposed  device  to  that  of  another  system 
that  observationally  has  appeared  to  be  stable. 

I  could  illustrate  the  mathematical  form  in  which  one 
might  pose  this  problem  in  our  application  by  considering 
explicitly  the  case  of  a  single  spatial  degree  of  freedom  -- 
such  as  the  transverse  particle  displacement, in  the  median 
plane  of  the  accelerator,  measured  from  the  equilibrium  orbit 

(Pig.  3). 

In  the  ideal  case  of  an  A-G  accelerator  that  is  perfectly 
constructed  and  aligned,  the  linearized  equations  of  motion 
would  be  a  simple  equation  of  the  form  of  Hill's  equation 
(to  which  the  Floquet  theory  would  apply)  --  as  shown  at  the 
top  of  the  Figure,  where  F(9)  is  periodic  in  9  and  in  practice 


might  be  represented  simply  by  one  of  the  two  forms  indicated. 


LINEARIZED  EQUATIONS 

^  ^  +  F(9)u  =  0 


d0 


where   F(9)  = 


'periodic  Square  Wave 
or 

a  +  b  cos  N9 


P 


TYPICAL  NON-LINEAR  EQUATIONS 
One  degree  of  Freedom 


^2  P  X  3 

.2_^  +  (a+b  cos  9)u  +  a(sin  N9)u''  -  P(cos  N9)u-^  =  0 


d9 


Two  coupled  degrees  of  Freedom 


2   2. 


^  +  (a+b  cos  N9)u  +  a(sin  NO) (u^-v^) -P(cos  N9)(u^-3uv  )=0 
d9 

2  2 

^  +  (a'-b  cos  N9)v  -  2a(sin  N9)uv  +  3P(cos  N9)u  v  =  0 


d9 


:2 


Fig.  3 


In  actuality,  some  non-linearity  of  the  restoring  (focusing) 
force  will  be  present,  either  deliberately  or  inadvertently 
introduced.   My  own  interest  in  the  effect  of  non-linearities 
arose  in  connection  with  designs  for  fixed-field  alternating- 
gradient  (FFAG)  accelerators,  as  proposed  by  Symon  and  others 


of  the  Midwestern  Universities  Research  Association  (M.U.R . A. ) j 
in  which  non-linearities  necessarily  would  be  prominent  at 
the  amplitudes  of  interest.   Representative  terms  that  could 
be  added  to  account  for  quadratic  and  cubic  components  to  the 
force  would,  in  that  application,  be  of  the  form  shown  in  the 
middle  of  the  Figure,  and  below  it  one  sees  the  Hamiltonian 
pair  of  equations  that  one  could  take  to  describe  the  corre- 
sponding motion  in  two  spatial  degrees  of  freedom. 

The  inclusion  of  a  periodic  9-dependent  coefficient  for 
the  linear  terms  in  these  equations  is  not  significant,  since 
that  feature  can  be  transformed  away  by  an  explicit,  well- 
behaved  transformation.   What  ^  important  is  that  we  do  not 
have  linear  equations,  to  which  the  well-behaved  properties 
of  the  Floquet  solutions  apply,  nor  do  we  have  a  9-independent 
Hamiltonian  that  could  be  taken  as  a  constant  of  the  motion. 
(As  is  of  course  known.  Dr.  Moser  has  investigated  the 
possibility  of  systematically  introducing  a  series  of  trans- 
formations for  the  working  variables  such  that  the  9-dependence 
of  the  Hamiltonian  function  becomes  displaced  to  higher  and 
higher  order  --  this  is  done  with  the  expectation,  as  I 
understand  it,  that  the  new  Hamiltonian  functions  will 
progressively  become  better  approximations  to  a  constant  of 
the  motion,  as  in  fact  certainly  seems  to  be  the  case  in 
limited  applications  of  this  technique  to  the  study  of  solutions 
over  moderately  long  intervals . ) 

One  has  the  impression  that  the  essential  features  of 
solutions  to  equations  of  this  form  can  be  exhibited  by  suitably 


constructed  non-linear  algebraic  transformations  whose 
short-term  phase-plane  characteristics  would  be  qualitatively 
similar  to  those  implied  by  the  differential  equations,  and 
for  which  computational  studies  would  be  both  faster  and  less 
affected  by  numerical  error.   Such  transformations  can,  in 
fact,  be  regarded  as  special  cases  of  the  differential 
equations,  with  the  periodic  coefficients  becoming  periodic 
delta  functions  of  the  independent  variable  (discrete, 
localized  lenses).   What  is  required  of  the  transformation  is 
simply  that  the  iterated  values  of  the  coordinate  and  conjugate- 
momentum  variables,  as  functions  of  the  previous  values  of 
these  variables,  satisfy  the  Poisson-bracket  conditions  -- 
or  any  other  of  the  several  alternative  forms  in  which  one 
can  express  the  necessary  and  sufficient  conditions  for  a 
system  to  be  Hamiltonian.   In  one  spatial  degree  of  freedom, 
this  condition  can  be  simply  stated  as  requiring  that  the 
transformation  be  area  preserving  in  the  two-dimensional   q,p 
phase  plane. 

A  simple  transformation  in  which  I  became  interested 
some  ten  years  ago  at  M.U.R.A.  pertained  to  motion  in  one 
spatial  degree  of  freedom  and  was  believed  to  simulate  the 
solution  to  the  radial  motion  in  the  median  plane  of  a  spirally- 
ridged  FFAG  accelerator  if  the  coordinate  and  momentum  for  this 
case  are  plotted  at  one-period  intervals  of  the  structure. 
This  transformation  is  shown  in  Fig.  4,  in  which  A  is  a 
constant  and  I  regard   x  and   y  respectively  as  a  coordinate 
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and  its  canonically-conjugate  momentum. 

X'    =  Ax  +    (l-A^)y  +    (1-A)[x  +    (l-A)y]^ 


y'    =   +  X  +  Ay   +   [x   +    (l-A)y] 


If     T(x,  ,y-,)    =    (X2,y2),      then     T   ^(x-^,-y^)    =    (x^^-yg) 


x-j-  ^  0.03061 


Sketch  for     A  = 


78 


cos   ^A  ^  (0.357^5)  {2Tr) 


Area  ^  5.5  x  lO""^ 
(would  shrink  to  zero 
if  A  -^  -  Vs) 


F.:  (-V:^5Iz5  ^  0)  =  (-0.053966.  ••,  0) 
^  26 

F^  .  (2-  ,   +  ^Ai-1^)  =  (0.03846---,  +  0.05688---) 
2'5   26   -   169 


Fig.  4 

The  upper  and  lower  signs  can  be  taken  as  referring  to  the 
forward  or  inverse  transformation,  respectively.   The  constant 
A  can  be  interpreted  as  the  cosine  of  the  phase  advance  per 
iteration  that  would  be  exhibited  by  solutions  to  the  linearized 
transformation,  and   A  normally  would  be  taken  to  have  an 
absolute  value  less  than  unity.   It  also  turns  out  that  to 


avoid  being  exactly  at  a  one-third  resonance  when  the  quadratic 
terms  are  present,  one  likewise  should  avoid  the  value 
A  =  -1/2   (for  which   cos"  A  =  27r/3 )  . 

This  transformation,  as  written,  was  intended  to  simulate 
the  solutions  in  a  spirally-ridged  FFAG  accelerator  when  viewed 
at  a  suitable  point  of  symmetry  in  the  structure  (and  at 
homologous  points,  spaced  at  intervals,  27r/N,  equal  to  one 
period  of  the  structure) .   In  this  form  the  transformation 
(T)   has  the  great  convenience  of  exhibiting  a  symmetry  about 
the  X-axis,  as  has  been  noted  on  the  Figure  and  as  is  evident 
from  inspection  of  the  transformation  equations.   The  phase 
diagram,  as  revealed  by  short-term  computations  --  or  roughly 
by  application  of  Moser's  methods  or  in  some  respects  by  more 
simple  analytic  considerations  --  looks  somewhat  as  sketched 
on 

and  in  this  respect  the  proposed  transformation  appears  con- 
sistent with  trajectories  computed  for  some  typical  FFAG 
structures  in  which  one  would  plan  to  employ  the  greater  part 
of  the  apparently  stable  region  that  is  situated  within  the 
roughly  triangular  boundary  shown. 

What  appears  to  be  a  stable  area  is  bounded  by  the  curves 
that  I  have  drawn  through  the  three  unstable  fixed  points  of 
order  3  (F-,  ,  Pp,  ^^  ^  '  that  are  rigorous  fixed  points  whose 
coordinates  can  be  explicitly  obtained  from  the  roots  of  a 
quadratic  equation.   In  calling  your  attention  to  transformations 
of  this  form  (and  to  the  analogous  non-linear  differential 


the  Figure  /for  A  =  -5/8,  cos  ^A  =  0.357^5  (2Tr)  =  128.?  deg_^7. 
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equations).   I  do  not  mean  to  imply  that  present  accelerator 
or  storage-ring  projects  visualize  such  pronounced  excursions 
toward  the  boundaries  of  regions  so  obviously  influenced  by 
non-linear  resonant  effects,  but  smaller  non-linearities  may 
well  be  present  in  practice  and  future  accelerator  concepts 
may  again  involve  strongly  non-linear  restoring  forces.   I 
also  would  not  wish  to  imply  that  the  most  important  problems 
to  pursue  in  this  connection  are  those  that  involve  only  one 
spatial  degree  of  freedom.   It  may  well  be  that  results  of 
some  computational  experiments  --  and  theoretical  work  by 
Kolmogorov  and  others  that  I  hope  Dr.  Moser  will  review  for 
us  --  can  set  the  designer's  mind  at  rest  with  respect  to 
possible  long-range  instability  in  one  degree  of  freedom  and 
that  emphasis  should  be  given  in  this  context  to  the  much 
more  difficult  problems  that  arise  with  systems  having  two 
(or  more)  spatial  degrees  of  freedom. 

In  any  event,  I  thought  last  spring  that  for  me  it 
would  prove  instructive  to  examine  the  performance  of  the 
transformation  just  shown,  using  the  CDC-660O  computer  in 
double  precision,  since  my  previous  look  at  this  problem  had 
been  with  the  earlier  IBM-704 .   This  work  was  undertaken  at 
Berkeley  with  the  programming  assistance  of  Mrs.  Levine  and 
with  helpful  advice  from  Eric  Beals  and  Loren  Meissner.   By 
using  double  precision,  we  had  available  some  96  bits  -- 
although  truncation  rather  than  true  rounding  was  used  in 
the  arithmetic  operations  --  and  it  was  possible  to  investigate 
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spurious  fluctuations  and  drifts  both  by  computational 
experiments  and/or  by  a  rough  analysis. 

In  investigating  the  behavior  of  an  algebraic  transfor- 
mation such  as  the  one  just  proposed,  one  of  course  has  a 
wide  choice  of  working  variables.   Thus,  if  desired,  new 
variables  could  be  introduced  that  would  result  in  a 
transformation  having  symmetry  about  the  +45   diagonal,  or 
one  could  select  variables  that  would  reduce  the  number  of 
computational  steps  (and  hence  improve  the  speed  and  accuracy 
of  the  work)  required  to  perform  each  iteration.   One  may 
merely  regard  the  introduction  of  such  new  variables  as  a 
mathematical  operation,  but,  in  some  instances,  use  of  the 
new  variables  rather  easily  can  be  interpreted  as  viewing 
the  dynamics  at  some  other  reference  point  within  the  lattice 
of  the  physical  structure. 

For  my  own  computational  experiments,  it  was  convenient 
to  introduce  new  ("working")  variables  defined  in  Fig.  5  f'^^ 
the  direct  or  inverse  transformations,  respectively.   In 
terms  of  these  variables  the  transformation  assumed  the  form 
shown  at  the  bottom  of  the  Figure. 
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By  introduction  of  the  working  variables 
q  =  X  +  (l-A)y 
P  =  +  y 

for  the  direct  transformation  T  or  its  inverse, 

q'  =  Cp  +  q 

p'  =  (q'-l)q'  +  p 
with 

C  =  2(1-A) 

Fig.  5 

This  transformation  seemed  particularly  suitable  for  repeated 
iteration,  but  phase  diagrams  plotted  in  terms  of  these 
working  variables  (q,p)  do  not  show  any  apparent  symmetry 
and  it  therefore  was  convenient  to  express  results  in  terms 
of  the  original  x,y  variables  whenever  output  data  were  printed 
The  parameter  A  (=  -5/8,  a  binary  fraction)  and  derived 
parameters  /such  as  (1-A)  =  13/8  and  C  =  2(1-A)  =  13/^7  were 
stored  exactly  in  the  computer  memory,  and  double-precision 
octal  I-O  was  available  if  required  for  resumption  of  any 
particular  run. 

An  initial  attempt  to  reveal  erratic  or  stochastic 
behavior  gave  negative  results,  in  that  for  runs  of  as  many 
as  10   forward  iterations,!^  starting  with  x  =  +0.026  =  O.85  x_, 
there  appeared  to  be  no  computationally  significant  anomolous 

]      9^1 -169  sec  of  CP  time  were  required  for  a  run  consisting 

7  7 

of  10   forward  iterations  followed  by  10   reverse  iterations 
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drift  --  e.g.,  after  correction  for  computational  drift  and 
allowance  for  computational  fluctuations,  the  +x  intercept 

of  the  apparent  invariant  phase  curves  could  be  said  to  remain 

-27 
constant  to  an  absolute  accuracy  of  about  10    ,  or  a  relative 

-26 
accuracy  of  about  4  x  10    in  this  example. 

A  positive,  and  hence  more  dramatic,  effect  can  be 
obtained,  however,  by  examining  in  some  detail  the  character 
of  the  apparent  separatrix.   For  this  purpose  we  note  that 
the  transformation  T    will  return  a  phase  point  to  any 
one  of  the  unstable  third-order  fixed  points  (F-,,  Fp,  F^ )  ^ 
and  that  the  " separatrices "  drawn  through  the  fixed  points 
can  represent  at  these  points  the  directions  of  outgoing  and 
ingoing  eigenvectors  for  the  transformation  T   when  the  latter 
is  linearized  about  the  fixed  points.   If  a  true  separatrix 
FpF^   exists,  a  line  segment  formed  of  points  situated  very 
near  to   Fp   and  that  lie  on  a  curve  whose  slope  is  that  of 
the  appropriate  (outgoing)  eigenvector  through  that  point 
should  iterate  under  T   so  as  to  approach  F^   along  a 
similar  eigenvector  direction  and,  in  the  process,  should 
generate  a  single  smooth  curve   (FgF-,).   The  symmetry  of  the 
present  transformation  (with  respect  to  the  x-axis)  would 
imply,  therefore,  that  the  curve  so  generated  cross  the  x-axis 
with  a  vertical  slope. 

Numerically,  forward  iteration  of  a  line  segment  from 
the  neighborhood  of  Fp   down  toward  the  x-axis  should  be 
favorable  with  respect  to  unavoidable  computational  errors. 


Ik 


and  detection  of  a  non-vertical  crossing  would  imply  the  non- 
existence of  a  firm  separatrix  between  Fp   and   F^   (and,  by 

2       -1 
application  of  T   and   T   or  T   ,   correspondingly  also 

would  imply  the  absence  of  firm  separatrices   F-zF-,   and   F.Fp.) 

Our  work  with  double-precision  arithmetic  made  it 
unambigously  clear  that  extension  of  the  eigenvector  direction 
at   Fo   led  to  non-vertical  crossing  of  the  x-axis  (at  crossing, 
slope  S  +10  ),   although  if  single-precision  arithmetic  had 
been  employed  the  results  might  have  been  ambiguous  in  this 
case.   /Checks  of  the  computational  process  could  be  made  by 
use  of  different  compiler  systems,  writing  the  algebraic 
statements  in  a  different  order  or  form,  influencing  the 
truncation  errors  (as  by  introducing  intermediate  steps  such 
as  "times  3"  followed  by  "divide  by  3"  into  one  or  more 
equations  for  each  iteration),  and  by  use  of  varied  starting 
conditions.   Although  some  of  these  steps  affected  noticeably  . 
the  location  of  individual  iterates  (x,y)  near  the  x-axis, 
they  had  no  significant  effect  on  the  curve  that  was  generated 
in  this  procesSj_7 

As  is  indicated  in  Fig.  6,  the  situation  thus  is  such 
that  a  line  segment  QRS  originating  at   Fp  with  the  eigenvector 
slope  will  transform,  after  repeated  iterations  of  T  ,   to 
generate  the  S-curve  shown  just  below  the  x-axis  (initially 
with  an  eigenvalue,  for   T  ,  A-,-  =  2.102). 
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y 


Fig.  6 

A  similar  sequence  of  inverse  iterations  from  F^   will 
generate  the  dashed  curve,  and  forward  iteration  wou3.d  cause 
points  on  the  dashed  curve  to  approach   F^   in  the  absence  of 
computational  errors.   Thus  only  particular  points  (e.g.: 
Q,  R,  and  their  iterates)  that  lie  on  both  curves  can  originate 
arbitrarily  close  to   Fp   and  ultimately  approach  F^ . .   The 
area  1  transforms  to  the  equal  area  3  (unit  Jacobian)  and 
areas  1  and  2  are  equal  by  virtue  of  the  symmetry  of  the 
transformation,  so  all  the  loops  shewn  --  and  their  iterates 
under  the  transformation  T  --  are  equal  in  this  example. 

It  will  be  recognized  that  points,  such  as   Q'   and   R', 
that  approach  F^   ultimately  do  so  in  smaller  and  smaller 
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steps   (^x-r  =  ''^T  ~  0.48).   Conservation  of  the  areas  mentioned 
thus  requires  that  the  curve  bounding  the  transforms  of  the 
areas  2  and  3  correspondingly  must  become  increasingly 
elongated  laterally,  and  (as  suggested  by  Fig.  7)  the  solid 
curve  Q'R'S'Q"  will  become  increasingly  sinuous. 


Fig.  7 

The  areas  shown  (transforms  of  regions  2  and  3)  each 

-11  8 

have  the  approximate  magnitude  4.8  x  10    ,  or  roughly  l/lO 

times  the  area  of  the  entire  superficially-stable  region 

FnFoF^.   The  uniqueness  of  the  transformation  precludes  that 

a  finger  such  as  that  shown  extending  to  the  left  will  be 

intersected  by  one  of  its  iterates,  although  such  an  iterate 

could  well  cross  the  dashed  line.   If  the  transformation 

continues  to  be  applied,  these  fingers  thus  cannot  remain 
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at  all  times  within  the  finite  area  that  at  first  might  have 
been  supposed  to  be  a  stable  region,  and  an  observed  means  of 
escape  is  that  in  which  an  iterate  of  this  finger  does  cross 
the  dashed  line  to  enter  a  loop  extending  to  the  right. 

It  will  be  appreciated  that  the  evolution  of  these 
fingers  necessarily  will  become  quite  complex,  and  certainly 
they  will  make  inroads  into  the  ostensible  "interior"  region, 
but  at  this  point  it  is  by  no  means  clear  that  the  entire 
"interior"  region  need  be  consumed  in  this  way.   What  is  clear 
is  that  the  non-smooth  joining  of  the  forward  and  backward 
eigenvector  directions  from  Fp   and  F^ ,   respectively  -- 
either  by  a  disparity  of  slope  (as  in  the  present  example) 
or  by  a  disparity  of  some  higher  derivative  --  wi.ll  imply  that 
these  curves  do  not  generate  a  true  separatrix  and  the  truly 
stable  area  (if  one  exists)  consequently  is  somewhat  smaller. 

In  the  "interior,  "  closer  to  the  "stable"  fixed  point 
at  the  origin,  other  (higher-order)  fixed-point  systems  may 
be  found.   It  can  become  increasingly  difficult  in  such  cases, 
however,  to  examine  whether  smooth  or  non-smooth  intersection 
of  their  extended  eigenvectors  occurs,  and  a  computational 
investigation  would  offer  little  hope  of  establishing  absolute 
smooth  joining. 

It  is  interesting  and  somewhat  informative  to  examine 
other  area-preserving  transformations  in  these  respects. 
Professor  deVogelaere,  who  is  now  at  Berkeley  and  has  been 
very  generous  in  discussion  of  these  points,  has  for  some  time 
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been  directing  attention  to  a  very  simple  transformation  that 
requires  little  calculational  precision  to  exhibit  the  failure 
of  eigenvectors,  when  extended,  to  join  smoothly.   In  its 
original  form  this  transformation  was  of  the  form  shown  at 
the  top  of  Fig.  8, 

De  Vogelaere's  Transformation 

2 
X '  =  y  +  X 

2 
y '  =  -X  +  X ' 

-^  Tr.  =  0. 

Modified  De  Vogelaere  Transformation 
X'  =  y  +  Tx  +  (l-T)x 
y'  =  -X  +  Tx'  +  (l-T)x'^ 
■^  Tr.  =  T. 

Fig.  8 

and  can  be  shown  to  possess  the  same  symmetry  property  with 
respect  to  the  x-axis  that  was  noted  for  the  transformation 
I  mentioned  originally.   This  particular  form  may  be  somewhat 
distressing  to  a  physicist,  since  the  trace  of  the  matrix 
for  the  linearized  transformation  vanishes  (phase  advance  =  7r/2), 
and  one  would  be  situated  at  a  "quarter-integral  resonance." 
A  simple  generalization  of  this  transformation  can  be  made, 
however,  without  affecting  the  interesting  features  by  adopting 
the  form  shown  at  the  bottom  of  the  Figure.   For  either 
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transformation,  an  unstable  fixed  point,  of  order  1,  occurs       ^ 
at  the  point  (1,0),  and  the  eigenvector  directions  extended 
from  this  point  show  a  gross  failure  to  join  smoothly  on 
crossing  the  negative  x  axis.   This  behavior   is  indicated 
on  Fig.  9. 

Although  this  transformation  may  not  be  representative 
in  any  evident  way  of  accelerator-orbit  behavior,  its  study 
can  be  informative  concerning  features  that  can  arise  in 
non-linear  problems.   One  can  find,  for  example  with 
T  =  -I/8,  systems  of  unstable  fixed  points  for  which  the  eigen- 
values are  negative  and  that  do  not  appear  to  have  associated 
with  them  a  corresponding  system  of  stable  fixed  points. 

/E.g.,  for  T  =  -1/8:   Order  4/1,  with  eigenvalues  -3.197 

and    -0.3128   at   X  =  O.5615983,    y   =  0; 

Order  8/2,  with  eigenvalues  -8.369 

and  -0.1195  at  X  =  0.4562733,  y  =  Oj 

Also,  with  T  >  0  (e.g.,  +1/8),  one  finds  for  example  a  system 

of  unstable  and  stable  fixed  points  of  order  19/4  in  the  "interior 

that  roughly  exhibits  the  classical  "island  structure"  but  for 

which  the  extended  eigenvectors  very  easily  can  be  shown  not 

to  intersect  smoothly  (both  on  the  inner  and  on  the  outer 

island  boundaries).   /One  member  of  the  stable  and  of  the         , 

—  1 

unstable  fixed-point  systems  of  this  family  will  be  found 

respectively  at  x  =  -0.38429776,  y  =  0  and  x  =  +0.50736937, 

y  =  0.   The  eigenvalues  for  the  latter  are  Aj  =  l/A-rj  ^3.221/ 
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U.S.F.P 


(1,0) 


H 


Fig.  9 

Somewhat  closer  to  the  origin,  however,  one  finds  a  similar 
system  of  order  9/2  (indicated  on  Fig.  lO) ,  for  which  the 
intersections  appeared  completely  smooth  to  the  limit  of  my 
computational  accuracy  /fixed-points  at  x  =  -0.29000009,  y  =  0 
and  at  X  =  +O.52176070,  y  =  0  (with  A^  =  iA-j-^  =  1.179)_7- 
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Fig.  10 

Experience  with  each  of  the  transformations  mentioned 
could  be  said  to  be  not  inconsistent  with  the  view  that  some 
of  the  curves  commonly  regarded  as  stability  boundaries  are 
not  firm  and  imply  a  sort  of  stochastic  behavior  of  phase  points 
in  their  neighborhood,  but,  so  far  as  we  can  tell,  others  may 
be  perfect  barriers  and  indeed  may  have  true  invarient  phase 
curves  nearby  (ignoring,  of  course,  the  fluctuations  of  the 
parameters  that  in  practice  would  be  present  as  a  result  of 
"noise"  in  a  physical  system).   The  designer  of  particle 
accelerators  --  or,  more  particularly,  of  storage  rings  -- 
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of  course  would  like  to  know  whether  any  regularity  of  behavior 
can  be  anticipated  in  this  regard,  and  whether  quantitative 
or  reasonably-accurate  semi-quantitative  estimates  can  be 
made  of  diffusion  rate  in  a  region  that  possibly  is  characterized 
by  stochastic  behavior.   It  may  be,  of  course,  that  such  detailed 
and  subtle  features  of  behavior  are  highly  sensitive  to  the 
exact  values  of  the  physical  parameters  that  determine  the 
transformation,  and  as  such  would  be  beyond  the  precise 
physical  control  of  the  experimenter. 

(I  might  add  that  the  systems  of  unstable  fixed  points 
that  have  been  found  for  both  my  transformation  and  for  that 
of  Professor   deVogelaere  appear  to  show  what  to  me  is  a  quite 
surprising  regularity:  I  believe  that  deVogelaere  noticed 
empirically  that  for  his  transformation,  linearized  about  the 
fixed  points  of  a  certain  class,  the  quantity  -^  Trace  -  1 
showed  an  exponential  dependence  on   the  fixed-point  order. 
It  may  be  that  other  interesting  features  of  such  transformations 
can  also  be  found  to  exhibit  regularities  --  such  as  that 
suggested  by  Mel'nikov  in  regard  to  the  dimensions  of  the 
areas  that  are  developed  by  the  failure  of  separatrices  to 
join  smoothly  --  and  that  some  of  these  regularities  might 
permit  one    to  draw  inferences  that  would  have  practical 
application. ) 

An  interesting  class  of  transformations  has  been  proposed 
by  Dr.  McMillan,  for  which  it  is  possible  to  establish  that  the 
eigenvectors  through  the  fixed  point  or  fixed  points  do  form 
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a  true  impenetrable  separatrix.   He  finds  it  convenient  in' 
t'.'iis  connection  to  consider  the  transformations  to  be  written 
in  a  form  showing  symmetry  about  the  445   axis  (rather  than 
about  the  x-axis)  --  a  form  in  which  each  of  the  transformations 
previously  mentioned  can  be  written  by  a  suitable  change  of 
variables.   This  form  is  taken  to  be,  for  the  forward  trans- 
formation, that  which  will  be  seen  on  Fig.  11  in  the  upper 
left-hand  corner.   Such  a  transformation  has  a  simple  physical 


X'  =  y 


y 


-X  +  f(y) 


Fig.  11 
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interpretation  in  terms  of  a  structure  formed  simply  from  a 
periodic  sequence  of  identical  non-linear  lenses   /Ay  =  -x  +  f(x_y7 
placed  with  unit  separation  just  before  the  measurement  planes  and 
with  linear  lenses  [Ey  =   -x/     situated  just  after  each  of  these 
planes.   For  a  fixed  point  of  order  one,  the  trace  of  the  matrix 
for  the  transformation  linearized  about  that  point  is  just  the 
derivative  of   f   at  that  point,  and  stability  of  the  linearized 


system  requires 


1  f, 


<  1. 

The  transformation  in  the  form  shown  can  be  executed  by  a 
simple  graphical  construction,  in  which  one  first  mirrors  the 
initial  point  about  the  diagonal  and  then  vertically  reflects 
the  resulting  y-value  about  the  curve   y  =  -^  f  (x)  .   A  fixed 
point  of  order  1  is  synonymous  with  the  intersection  of  the 
diagonal  with  this  latter  curve.   Referring  to  Fig.  12, 
McMillan  has  pointed  out  /UCRL-177957  that  if  the  function 
f(y)  can  be  expressed  as   0(y)  +  0~  (y),   where   0~   is  the 
inverse  of  0,   then  x  =  0(y)   will  be  an  invariant  curve. 
The  phase  diagram  shown  in  the  Figure  is  for  an  example  in  which 


X  =  0(y)  =  a ^  and    f (y)  =  g  ^^  g   , 

y  +  a  a   -  y 

with  fixed  points  at   x  =  y  =  +  y/a   -  1    (|a|  >  l) . 
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Example  in  which 

f(y)  =  0(y)  +  0~^{y) 


for   a  =  -^/4] 


Fig.  12 

An  interesting  and  simple,  although  perhaps  artificial, 
example  of  this  class  is  one  suggested  by  Dr.  Judd,  for  which, 
as  indicated  on  Fig.  13 >  one  employs  a  step-wise  linear  f (x) . 
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Fig.  13 


Here  the  "broken-line  separatrlx  shown  is  a  rigorous  barrier. 
Also,  of  course,  motion  of  sufficiently  small  amplitude  to 
remain  within  the  dotted  square  will  trace  the  perfect 
elliptical  invariant  phase  trajectories  that  are  characteristic 
of  truly  linear  motion.   Outside  of  this  dotted  square,  moreover, 
stable  fixed  points  can  be  found   /for  example,  with  a  =  0.6: 
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an  order-6  FP  at  x  =  y  =  9/13  and  an  order-7  FP  at 
-X  =  y  =  378/6837  that  are  surrounded  by  an  elliptical 
boundary  within  which  the  motion  is  again  truly  linear  (and 
truly  stable  "in  the  small'*).   In  Fig.  l4  one  such  region  is 
shown  as  the  hatched  area  in  the  upper  right-hand  portion  of 
the  diagram  --  the  motion  being  controlled  in  such  a  case  in 
a  periodic  way  by  definite  values  of   f '   at  successive 


STABLE  AREA 
with  fixed  point 
of  Order  6- 


Fig.  1^ 
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iterations.  Other  points,  including  those  neighboring  unstable 
fixed  points  of  the  same  family,  however,  show  the  failure  of 
extended  eigenvector  directions  to  intersect  smoothly  and 
evidently  are  characterized  by  a  sort  of  stochastic  behavior 
that  is  confined  to  a  finite  region. T  Nonetheless,  the  motion 
within  the  diamond  is  not  completely  ergotic,  since,  as  we 
have  seen,  there  exist  in  this  area  regions  of  limited  extent 
within  which  the  motion  follows  true  simple  invariant  phase 
trajectories . 

In  presenting  this  material  I  do  not  mean  to  imply  that 
it  contains  much  information  that  is  new  to  the  mathematicians 
who  have  been  following  recent  developments  of  the  theory. 
I  am  afraid  that  some  of  the  computational  results  fail,  as 
might  be  expected,  to  give  definitive  experimental  answers 
to  some  of  the  significant  questions,  and  my  own  work  has 
neglected  entirely  the  more  important  (and  more  difficult) 
problems  that  arise  with  motion  in  more  than  one  spatial 
degree  of  freedom.   I  also  would  not  wish  to  claim  that  the 
questions  of  long-range  stability  are  of  quite  as  immediate 
importance  to  the  accelerator  designer  as  those,  for  example, 

^   In  this  transformation  (and  also  in  a  similar  trans- 
formation with   f(y)   analytic),  some  cases  were  found  in 
which  the  eigenvectors,  when  extended,  generated  a  double-S 
(or  double-Z),  rather  than  a  single  S  (or  Z)  between 
successive  iterates. 
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that  concern  collective  motion.   I  would  like  to  suggest, 
however,  that  the  properties  of  such  transformations  as  the 
one  I  mentioned  at  the  beginning  suggest  questions  that 
inherently  are  quite  fascinating,  and  that  it  would  be  quite 
helpful  for  future  design  work  if  our  mathematical  colleagues 
could  assist  us  in  obtaining  a  better  --  and  hopefully  somewhat 
quantitative  --  understanding  of  these  properties. 
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Appendix  A 

I  would  like  to  call  specifically  to  your  attention  a 
situation  we  encountered  in  which  an  area-preserving  algebraic 
transformation,  when  linearized  about  certain  high-order  fixed 
points,  exhibited  negative  eigenvalues.   This  phenomenon  arose 
in  connection  with  some  runs  that  I  made  using  a  modification 
of  Professor  deVogelaere ' s  transformation,  namely 

2 
/'  X'  =  y  +  Tx  +  Cx 


T; 


2 
y'  =  -X  +  Tx'  +  Cx' 


with   C  =  1  -  T   and   T  =  -1/8. 

This  transformation,  whose  phase  plots  are  symmetrical 
about  the  x-axis,  has  an  unstable  fixed  point  of  order  1  at 
the  point  (1,0)  and  the  eigenvector-directions  extended  from 
this  point  develop  loops  in  the  familiar  way  (as  qualitatively 
sketched  in  Fig.  1).   In  seeking  higher-order  fixed  points  in 
the  "interior"  of  this  diagram,  the  following  families  of 
unstable  fixed  points  appeared  for  which  the  eigenvalues  are 
negative : 

Family  1  ("Tune"  =  1/4): 

X  y 


0.5615983  0 

0.2846170  -0.5060428 

-0.4504872  0 

0.2846170  0.5060428 

Eigenvalues:    A   =   -3-197,    -0.3128 
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Fig.    1 
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Family   2    ("Tune"=   2/8): 

X  y 


0.4562733  0 

0.1771744  -0.4431055 

-0.4299377  0.0845201 

0.3462146  0.5215085 

0.6130793  0 

0.3462146  -0.5215085 

-0.4299377  -0.0845201 

0.1771744  0.4431055 

Eigenvalues:  A  =  -8.369,  -0.1195 

Phase  trajectories  formed  by  extending  the  eigenvector 
directions  for  the  system  described  as  Family  2  are  shown  in 
Fig.  2  in  the  region  near  the  points 

(0.6130793,  0)   and   (0.4562733,  0). 

It  does  not  appear  necessary  to  have  an  associated 
family  of  stable  fixed  points  in  such  a  negative-A  case,  and 
none  was  found . 

With  reference  to  Fig.  2,  we  may  consider  a  point  such 
as   A  ,   that  lies  at  the  intersection  of  a  phase  trajectory 
that  approaches  the  right-hand  one  of  the  two  fixed  points 
shown  and  a  trajectory  that  is  directed  away  from  the  left- 
hand  point  in  the  figure.   After  8  applications  of  the  trans- 
formation, this  point  moves  so  as  to  lie  closer  to  the  right- 
hand  fixed  point  (but  on  the  extended  eigenvector  segment  that 
lies  on  the  opposite  side  of  this  fixed  point,  since 
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Fig.    2 
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A  =  -O.II95  "^  0),   while  it  also  comes  to  lie  more  remotely 
from  the  left-hand  fixed  point  (and  on  the  extension  of  the 
eigenvector  segment  on  the  opposite  side  of  this  fixed  point, 
since   A  =  -8.369  *=   O)  .   This  new  point  is  designated   A   0 
on  the  drawing. 

The  point   B    lies  on  the  intersection  of  outward-  and 
inward-directed  eigenvectors  from  the  left-hand  fixed  point. 

o 

Its   T    iterate  is  seen  to  lie  on  another  such  intersection  -- 

closer  in  on   the  eigenvector  for  which   |a|  <  1   and  more 

remotely  on  the  extension  of  the  vector  for  which   |a|  >  1. 

This  performance  may  suggest  that  it  is  not  necessary  for  a 

loop  such  as  that  shown  between  B   and  the  left-hand  fixed 
^  n 

point  to  enclose  another  (e.g.,  stable)  fixed  point  of  this 
rotation  number.   I  thus  wonder  whether  the  behavior  of  such 
fixed  points  (with   A   real  and  negative)  may  be  qualitatively 
distinct  from  that  expected  in  cases  for  which   A  >  0   (and 
for  which  one  expects  to  find  a  "string  of  islands" such  that 
one  then  inquires  concerning  the  smoothness  or  lack  of  smoothness 
with  which  eigenvector  directions,  when  extended,  intersect). 
The  type  of  phase  trajectories  indicated  in  Fig.  2  may  be  most 
readily  adaptable  to  fit  in  with  the  rather  gross  loops  that 
develop  in  this  region  of  the  phase  plot  from  the  eigenvectors 
extended  from  the  order-1  fixed  point  at  (1,0).   [To  avoid 
confusion  on  the  drawing  of  Fig.  2,  I  have  not  shown  there 
the  order-'^-  unstable  fixed  point   (O.5615983,  O),  of  Family  1, 
that  falls  in  this  region  of  the  diagram.   This  point  does  not 
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lie  on  the  intersection  of  any  of  the  Family-2  trajectories 
shown  in  the  Figure  --  it  lies  on  the  x-axis  some  three-and- 
one-half  or  four  millimeters  to  the  right  of  the  point   A  . 
The  features  illustrated  in  this  diagram  can  be  exhibited 
without  any  exceptional  computational  accuracy,  and  I  have 
found  it  convenient  to  ad  lib  runs  by  sitting  at  a  readily- 
programmed  Olivetti-Underwood  "101  Programma".] 

As  a  second  point,  I  neglected  to  mention  on  Tuesday, 
during  the  discussion  concerning  "my"  transformation,  that 
this  transformation  also  has  an  order-1  fixed  point  at  (1,0). 
The  transformation  in  question  had  the  form 

X'  =  Ax  +  (1  -  A^)y  +  (1  -  A)[x  +  (1  -  A)y]^ 

2 

y'  =  -x  +  Ay  +  [x  +  (1  -  A)y]  . 

The  extensions  of  the  eigenvectors  from  this  order-1  fixed 
point  also  develop  gross  loops,  as  might  be  expected.   This 
behavior  is  sketched  in  Fig.  3^  where  the  hatched  area  is  the 
triangular  area  (with  order-3  fixed  points  at  its  vertices) 
to  which  I  directed  attention  in  my  talk.   As  before,  the 
constant   A  has  been  given  the  value   -5/8. 

This  same  transformation,  but  with  A  =  -l/^,  shows  a 
quite  striking  approximate  constancy  of  the  "tune"  with  respect 
to  amplitude.   With  this  value  of  A  (A  =  -l/^)   the  small- 
amplitude  tune  corresponds  to  about  0.2902  {2jr)      radians 
phase  advance  per  iteration,  and  a,  system  of  7  stable  and 
7  unstable  fixed  points  (tune  =  2/7  =  0.2857)   develops  at  a 
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A 


-5/8 


Fig.    3 
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rather  large  amplitude.   [The  x-coordinates  of  the  fixed 
points  that  lie  on  the  x-axis  are  given  below: 

Stable  F.P.:   -0.123211  58766  12485  36429  04o273; 

U.S.  F.P.:     +0.2  (exact),  (half  trace  +  I.05078125). 

Reasonably-appearing  phase  curves  can  be  drawn  inside  and 

outside  tliis  system  (Fig.  4).] 

As  a  final  point,  you  may  recall  tha,t  during  my  talk  I 

mentioned  in  passing  that  certain  regularities  appeared  to  be 

present  in  the  eigenvalues  (or,  more  simply,  in  the  matrix 

trace)  for  the  successive  families  of  fixed  points  of  a 

particular  class.   You  might  be  interested  in  the  results  for 

fixed  points  of  tune  — for  Professor  deVogelaere ' s 

4m  +  1 
transformation  (the  transformation   T,   with  the  parameter 

T   set  equal  to  zero).   These  results,  from  the  CDC-660O, 

are  given  in  the  attached  Table,  titled  "TRACE",  and  apply 

to  this  transformation  after  linearization  about  the  fixed 

points  in  question. 

It  will  be  noted  that  the  amount  by  which  -^  Trace 

differs  from  unity  is  closely  the  same  in  absolute  value  for 

the  stable  and  unstable  fixed-point  families  of  the  same  order, 

particularly  for  those  with  a  tune  having  a  large  denominator. 

As  Professor  deVogelaere  has  pointed  out,  a  plot  of 

log|l  -  i  Trace]  vs  this  denominator  is  highly  linear  through 

many  decades.   Somewhat  similar  regularities  also  appeared 

with  my  transformation  (taking  A  =  -5/8) . 
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DISCUSSION  FOLLOWING  LASLETT'S  PAPER 


Jurgen  Moser* 


I  should  like  to  add  some  comments  about  some  recent 
theoretical  work  on  the  problems  of  the  type  discussed  by 
Laslett  as  well  as  mention  some  numerical  experiments 
conducted  by  other  people. 

It  is  hard  to  believe  that  these  questions  have  their 
origin  in  the  theory  of  accelerator  design.   In  fact,  orig- 
inally the  basic  principle  for  the  orbit  stability  of  the 
Alternating  Gradient  Synchrotron  depends  to  a  large  extent  on 
the  linear  theory.   However,  in  the  course  of  deeper  studies 
of  the  particle  orbits  in  an  accelerator,  the  nonlinear 
effects  became  more  and  more  important  and  a  source  of  concern. 
It  was  crucial  to  have  some  guarantee  that  the  nonlinear  terms 
would  not  destroy  the  stability  which  was  so  carefully  provided 
for  by  the  linear  theory.   This  led  to  very  delicate  and 
difficult  theoretical  problems  which  also  were  approached  with 
numerical  experiments  conducted  on  simplified  model  equations 
as  we  have  just  learned  from  Dr.  Laslett 's  presentation. 

First,  to  mention  some  numerical  studies  of  quadratic 

measure  preserving  mappings,  I  want  to  refer  to  the  interesting 

work  of  Henon,  "Numerical  Study  of  Quadratic  Area-Preserving 

Mappings",  1967,  to  be  published.   He  studied  simple  quadratic 

mappings  of  the  same  nature  as  Laslett  mentioned  toward  the 

end  of  his  talk  (transformation  by  Dr.  McMillan).   The  main 

*  Courant  Institute  of  Mathematical  Sciences,  New  York 
University 
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results  are  the  following: 

In  the  neighborhood  of  an  elliptic  fixed  point,  one 
finds  well  organized  curve  patterns  which  disintegrate  at 
larger  distances.   Immediately  the  problem  arises  to  determine 
the  size  of  this  "stability  region"  as  well  as  determine  this 
curve  pattern  possibly  by  analytic  procedure.   Henon  restricted 
himself  to  the  description  of  the  computations  and  comparison 
of  his  results  with  asymptotic  series  obtained  by  G.  D.  Birkhoff, 
The  remarkable  fact  is  that  the  agreement  between  the  numerical 
results  and  these  asymptotic  series  is  tremendously  close. 
However,  no  clear  indication  for  the  size  of  this  stability 
region  is  in  sight. 

It  would  require  more  space  to  refer  to  various  theoret- 
ical studies  which  have  been  made  in  this  direction.   First 
of  all,  it  is  to  be  mentioned  that  the  curve  pattern  observed 
in  most  of  these  problems  in  general  does  not  exist  in  the 
mathematically  rigorous  sense  as  a  continuous  family  of 
invariant  curves.   This  fact  was  known  to  Poincare  and  is 
closely  related  to  his  nonexistence  proof  of  integrals.   Of 
course,  one  can  easily  construct  examples  for  which  such 
family  of  closed  curves  does  exist;  however,  this  is  an 
exceptional  situation.   For  a  rigorous  proof  of  this  last 

statement  we  refer  to  Rllssmann  .   The  situation  is  actually 
^ 

1   H.  RUssmann,  'HJber  die  Existenz  einer  Normalform  inhaltstreuer 
elliptischer  Transf ormationen,"  Math.  Ann.  137,  196?, 
pp.  64-6?. 


50 


very  delicate  insofar  as  arbitrarily  small  changes  of  the 
mapping  lead  to  a  distribution  of  the  curve  pattern.   On 
the  other  hand,  by  changes  of  sufficiently  high  order  terms, 
one  can  force  the  existence  of  a  family  of  smooth  closed 
curves.   The  situation  is  somewhat  reminiscent  of  the  distri- 
bution of  rational  numbers  on  the  real  line.   Both  the 
rational  numbers  on  the  real  line.   Both  the  rationals  as 
well  as  the  irrationals  are  dense  yet  the  irrationals  form 
the  majority. 

To  continue  this  analogy:   it  is  rather  easy  to  establish 
the  existence  of  irrational  numbers  and,  similarly,  it  is  not 
too  hard  to  establish  the  existence  of  mappings  which  do  not 
possess  such  a  continuous  curve  pattern.   It  is  usually  a 
more  difficult  problem  to  decide  whether  a  preassigned  number, 
say  TT   is  irrational.   Therefore,  one  would  expect  it  should 
be  difficult  to  decide  whether  a  given  mapping  of  this  sort 

belongs  to  one  class  or  the  other.   In  this  connection,  I 

2 
want  to  refer  to  a  paper  of  mine   in  which  it  is  shown 

rigorously  that  even  the  simple  polynomial  mappings  do  not 

possess  a  family  of  closed  curves.   This  includes,  in 

particular,  the  example  by  McMillan. 


2 


J.  Moser,  "On  the  Integrability  of  Area-preserving  Cremona 
Mappings  near  an  Elliptic  Fixed  Point,"  Bol.  Soc.  Mat. 
Mexicana,  1961,  pp.  176-I80. 
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At  this  point,  it  seems  confusing  that  the  numerical 
evidence  does  not  match  the  theoretical  prediction  of  non- 
existence of  such  a  curve  pattern.   It  turns  out  that  this 
paradox  can  be  resolved:   even  though  there  is  no  continuous 
family  of  closed  invariant  curves  there  does  exist  a  large 
set  of  closed  curves  which,  however,  do  not  join  together 
into  a  continuous  family.   They  are  interrupted  by 
infinitely  many  gaps  so  that  one  obtains  a  "Cantor  set"  of 
closed  curves.   This  was  proven  in  my  paper   as  well  as  in 
Arnold's  paper  .   It  is  to  be  emphasized  that  the  closed 
curves  so  obtained  are  by  no  means  pathological,  they  are 
differentiable,  and,  in  fact,  analytic  curves.   However, 
their  distribution  is  pathological.   The  crucial  point  of 
these  investigations  is  that  they  insure  stability  of  the 
elliptic  fixed  point  since  points  starting  in  the  interior 
of  such  curves  are  forced  to  remain  there. 

What  happens  in  the  gaps  left  by  the  invariant  curves 
mentioned  above  is  only  partially  understood.   These  regions 
contain,  in  general,  infinitely  many  fixed  points  of  the 

■^      J.  Moser,  "On  Invariant  Curves  of  Area-preserving  Mappings 
of  an  Annulus, "  Nachr.  Akad .  Wiss.,  CWttingen,  Math.  Phys. 

Kl.  Ila,  Nr.  1,  1962,  pp.  1-20. 

2i 

V.  I.  Arnold,  "Proof  of  A.  N.  Kolmogorov ' s  Theorem  on  the 

Preservation  of  Quasi-periodic  Motions  under  Small 

Perturbations  of  the  Hamiltonian ",  Uspekhi  Mat.  Nauk 

U.S.S.R.  18,  Ser.  5  (113),  1963,  pp.  13-^0. 
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mapping  or  one  of  its  higher  iterates.   However,  the  behavior 
of  the  consecutive  iterates  in  these  regions  is  quite  irratlc 
and  irregular,  in  particular  near  unstable  fixed  points.   These 
domains  are  often  referred  to  as  "regions  of  instability".  For 
an  instructive  schematic  picture  we  refer  to  the  paper  by 
Arnold  .   There  are  conjectures,  more  or  less  informal,  that 
these  regions  of  instability  contain  open  invariant  sets  in 
which  the  mappings  are  ergodic .   Nothing  of  this  sort  has 
been  proven  and  there  is  some  numerical  evidence,  especially 
in  the  work  of  John  M.  Greene,  Forrestal  Research  Center, 
Princeton,   for  extremely  ergodic  behavior  in  these  domains. 
In  fact,  he  expressed  the  opinion  that  there  are  invariant 
domains  in  which  the  unstable  fixed  points  are  dense.   Since 
this  situation  is  typical  for  ergodicity,  one  may  take  this 
as  an  indication  for  the  existence  of  ergodic  domains. 

Comparing  the  numerical  and  theoretical  results,  one 
has  to  be  very  careful:   the  theoretical  results  are  strictly 
applicable  only  in  extremely  small  neighborhoods  of  elliptic 
fixed  points  (in  some  paper,  Henon  estimates  a  radius  of 
validity  of  the  proof  as  about  lO'^'^'^l).   This  is  due  to  the 


^   V.  I.  Arnold,  "Small  Divisor  and  Stability  Problems  in 
Classical  and  Celestial  Mechanics, "Uspekhi  Mat.  Nauk,  l8. 
No.  6  (114),  1963,  pp.  81-192. 

Private  communication. 
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shortcoming  of  proofs  which  are  very  crude  in  the  details  of 
the  estimate.   One  would  have  to  resort  to  numerical  experiments, 
such  as  Laslett  presented  here,  to  obtain  the  feeling  for  the 
actual  size  of  the  stability  region.   In  all  these  numerical 
studies  one  finds, in  fact,  a  sizeable  neighborhood  covered 
to  a  large  extent  by  invariant  curves.   On  which  parameters 
or  properties  the  size  of  this  stability  region  depends  seems 
still  to  be  quite  in  the  dark.   For  this  reason,  I  consider 
the  closer  study  of  numerical  experiments  extremely  worthwhile 
and  instructive. 
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ADDENDUM  TO  UCRL  17795,  "SOME  THOUGHTS  ON  STABILITY  IN 
NONLINEAR  PERIODIC  FOCUSING  SYSTEMS",  Sept.  5,  1967- 

Edwin  M.  McMillan* 

March  29,  1968 


1.)   Introduction 

In  UCRL  17795,  it  was  shown  that  curves  in  the  x,y 
plane  having  reflection  symmetry  about  the  positive  diagonal 
are  invariant  under  the  transformation: 

X'  =  y 

(1). 

y'  =  -x+f(y), 

where   f(y)   is  the  sum  of  the  two  values  of  x  corresponding 
to  the  given  y.   It  is  required  that  there  be  just  two  values, 
but  the  two  branches  on  which  they  occur  are  not  required  to 
have  a  common  analytic  form.   An  example  given  was  the  pair 
of  rectangular  hyperbolas   y  =  1  -a/(x+l)   and   y  =  -1  +a/(l-x), 
with   f(y)  =  2  ay/(l-y  ),   mentioned  in  paragraph  3  and  illus- 
trated in  Fig.  1.   The  question  whether  there  are  other 
invariant  curves  belonging  to  the  same   f(y)   was  left  open. 

This  question  was  answered  by  John  M.  Greene  in  a  letter 

to  L.  Jackson  Laslett  (March  8,  1968) .   He  pointed  out  that 

2      2 
all  curves  of  the  form   (1-x  )(l-y  )  +  2  axy  =  const.   are 

2 
such  invariants.   If  the  constant  has  the  value   2a-a  ,  the 

equation  factors  into  two  equations  representing  the  rectangular 

hyperbolas,  which  are  now  seen  to  be  simply  the  separatrices 

Lawrence  Radiation  Laboratory,  University  of  California, 
Berkeley,  California 
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of  a  family  of  invariant  curves.   In  the  course  of  checking 
the  invariance  of  "Green's  function"  by  the  methods  of 
UCRL  17795,  I  found  that  it  is  a  special  case  of  a  broader 
class,  which  can  be  called  "double  quadratic"  curves. 
2.)  "Double  quadratic"  curves 

Any  equation  which  is  quadratic  in  x  can  be  solved 
explicitly  for   x.   If  x  and  y  occur  in  it  symmetrically, 
it  represents  a  curve  with  the  required  symmetry  about  the 
positive  diagonal.   The  most  general  equation  with  these 
properties  is: 

Ax^y^  +  B(x^y+xy^)  +  C(x^+y^)  +  Dxy  +  E(x+y)  +  F  =  0,     (2) 
whose  solution  is: 


x  = 


2(Ay^  +  By  +  C) 


(By^+Dy+E) 


O) 


+ 


y(By^+Dy+E ) ^ -4 ( Ay ^+By+C ) ( Cy ^+Ey+F ) 


The  sum  of  the  two  values  of  x  gives   f(y) 


f(y)  .  _  By^  +  Dy  +  E  (2^) 

Ay^  +  By  +  C 


Since   f(y)   does  not  depend  on  F,   all  members  of  the  family 
generated  by  giving  different  values  to  F  are  invariant  under 
the  transformation  (l),  with  f(y)   given  by  (4). 
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We  thus  have  the  remarkable  result  that  an  f(y)   which 
is  the  ratio  of  any  two  quadratic  functions  of  y  leads  to 

a  family  of  invariant  curves,  with  the  single  restriction 

2 
that  the  coefficients  of  y    in  the  numerator  and  of  y   in 

the  denominator  must  be  of  equal  magnitude  and  opposite  in 

sign. 

The  first  order  fixed  points,  if  they  exist,  are  at 

f(y)  =  2y,   and  are  therefore  the  solutions  of: 

2  Ay^  +  3  By^  +  (2C+D)y  +  E  =  0  (5). 

The  number  of  parameters  in  (4)  is  easily  reduced; 
E  can  be  eliminated  by  a  coordinate  displacement  along  the 
positive  diagonal,  either  A   or  B   can  be  made  equal  to 
D   or   E  by  a  change  of  scale,  and  any  one  of  the  remaining 
parameters  can  be  set  equal  to  unity.   Thus  we  have  a  two- 
parameter  system.   Some  interesting  cases  are: 

(1)      A  =  1,  B  =  0,  C  =  -1,  D  =  2a,  E  =  0,  F  =  c. 

2    2  2      2 

xy     -x-y     +2axy  +  c=0.  ("Green's   function") 

f  (y)    =   ^^2 

1  -  y 

The  first  order  fixed  points  are  at  y  =  0,  +  \/l-a. 

The  separatrices  are  displaced  rectangular  hyperbolas,  as 

pointed  out  above. 
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(2)  A  =  1,  B  =  0,  C  =  1,  D  =  -2a,  E  =  0,  F  =  c. 

2  2     2   2 
xy  +x+y  +2axy+c=0. 

f  (y)  =  i-±^ 

1  +  y 

The  first  order  fixed  points  are  at   y  =  o,  +  j/a-1 

The  separatrix  is  the  curve  given  by  setting   c  =  0. 

In  cases  (l)  and  (2),  if  a   is  negative,  the  curve  is  rotated 

by  90°,   and  the  first  order  fixed  points  (except  the  one 

at   X  =  O)   become  second  order  fixed  points.   (See  paragraph  6 

and  Fig.  3b  of  UCRL  17795) 

(3)  A  =  0,  B  =  1,  C  =  -1,  D  =  0,  E  =  0,  F  =  c. 

2       2     2     2 
xy+xy   -X   -y  +c=0. 

f(y)  -  ^  - 


Th 


1  -  y 

2 


e  first  order  fixed  points  are  at   y  =  0,  ^ 

o 

The  separatrices  are  the  curve  given  by  setting   c  =  —  , 

27 
the  line  x+y+2  =  0,   and  the  curve   xy  -x-y  +2=0. 

(I  thank  Dr.  Laslett  for  finding  the  last  two  of  these.) 

(4)      A  =  1,  B  =  -2,  C  =  1,  D  =  0,  E  =  0,  F  =  c. 

P  p      p     p     P   ? 
X  y^  -2  (x^y+xy  )  +  x  +y   +  c  =  0. 

f(y)--^ 
(1-y)^ 

The  first  order  fixed  points  are  at   y  =  0,  p  (3  +  \/5)  • 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  Implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  Information  contained  In  this  report, 
or  that  the  use  of  any  Information, 
apparatus,  method,  or  process  disclosed 
In  this  report  may  not  Infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  Information,  apparatus, 
method,  or  process  disclosed  In  this 
report. 

As  used  In  the  above,  "person  acting  on  behalf 
of  the  Commission"  Includes  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  infor- 
mation pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with 
such  contractor. 
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